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Abstract 

I A state-dependent 1-limited polling model with A'^ queues is analyzed. 

. The routing strategy generalizes the classical Markovian polling model, in 

the sense that two routing matrices are involved, the choice being made 
according to the state of the last visited queue. The stationary distribution 
■ of the position of the server is given. Ergodicity conditions are obtained 

by means of an associated dynamical system. Under rotational symmetry 
assumptions, average queue length and mean waiting times are computed. 

> ! 1 Introduction 
m 

I Consider a taxicab in a city in which there are N stations at which clients arrive 

• and wait for the vehicle. When their turn has come to be served, they ask for a 



o 



transit to a destination (one of the other stations) where they leave the system. 
Whenever the taxicab finds a station empty, it goes somewhere else to look for 
I a client. The choice of the destinations by a client or by the empty taxicab are 

made via the two distinct routing matrices P and P. 

This system can also be seen as a polling model with N queues at which cus- 
tomers arrive and one server which visits the queues according to the following 
I rules: if the server finds a client at the current queue, it serves this client and 

chooses a new queue according to the routing matrix P; otherwise it selects the 
next queue according to P. This polling scheme is an extension of the classical 
Markovian polling model, with routing probabilities depending on the state of 
the last visited queue. As a taxicab only takes one client at a time, the service 
strategy is the 1-limited strategy, which is known to be more difficult to analyze 
than either the gated or exhaustive strategies — in which the server respectively 
serves the clients present on its arrival or serves clients until the queue is empty. 

Since the bibliography on polling models is plethoric, we refer the reader to 
the references given in Takagi [16] and Levy and Sidi [10]. Among the studies 
related to our work, we can cite Kleinrock and Levy [9], who compute waiting 
times in random polling models in which the destination of the server is chosen 
independently of its provenance, and Boxma and Weststrate [2] who give a 
pseudo-conservation law for a model with Markovian routing. Ferguson [6] and 
Bradlow and Byrd [3] study approximatively a model where the switching time 
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is station-dependent. Srinivasan [15] analyzes a polling system with the same 
routing policy as in our model, and with different service policies at each queue. 
However, the waiting times are not computed in the case where the routing is 
state dependent. Another important issue concerns state-independent ergodicity 
conditions. In this case, necessary and sufficient conditions were obtained by 
Fricker and Jai'bi [7, 8] for deterministic and Markovian routing with various 
service disciplines. Borovkov and Schassbergcr [1] and FayoUc et al. [5] give 
necessary and sufficient conditions for a system with Markovian routing and 
1-limited service. Finally, Schassberger [13] gives a necessary condition for the 
ergodicity of a polling model with 1-limited service and a routing that depends 
on the whole state of the system. 

Here, we give a new method to get the stationary distribution of the position 
of the server and, in the case of a fully symmetrical system, a way to compute 
the mean waiting time of a customer. This extends known results on 1-limited 
polling models, especially for the symmetrical case. However, we do not provide 
a pseudo-conservation law as for example Boxma and Wcststratc [2], since the 
computations are not as simple as in the state-independent routing case. Our 
method of proof allows to make only minimal assumptions on the arrival process 
and applies to cither discrete or continuous time models. Moreover, in the 
symmetrical case the results allow to compare various polling strategies. 

The paper is organized as follows: in Section 2, we present the model and 
give a functional equation describing its evolution. In Section 3, stationary 
probabilities of the position of the server are obtained and Section 4 is devoted 
to the general problem of ergodicity. Section 5 presents formulas to compute the 
first moment of the queue length at polling instants. These results arc used in 
Section 6 to calculate the waiting time of an arbitrary customer. Applications 
to known models are also given. 

2 Description of the model 

The system consists of N stations at which clients arrive according to a station- 
ary process. Let <S = {1, . . . , A''} be the set of stations. Assume that the server 
arrives after n — 1 moves at station i € S where a client is waiting. Then the 
server loads this client and goes to station j G S with probability pij, such that 
Pi,i + • • • + Pi,N — 1- Conversely, if station i is empty, the server polls station 
j with probability pij. The service policy is known as 1-limited policy, since 
at most 1 customer is served each time a station is visited. The two transition 
matrices will be denoted respectively by 

P = {Pi,j)i,j€S, P = {Pi,j)i,jeS- 

The number of new customers arriving to station q G S between the polling 
of stations i and j when there have been a service (resp. no service) at station 
i is Bi^.q{n) (resp. Bij.q{n)). The vectors Bi^j{n) = {B.,^j.i{n), . . . , Bij.N{n)) 
(resp. Bij{n)) are i.i.d. for different n, but their components may be dependents 
and not identically distributed. In the classical polling terminology, Bij{n) is 
the number of arrivals during both the service and the switchover and Bij(n) is 
the number of arrivals during a switchover. However, our model covers a wider 
class of applications. In particular, the service time may depend on both i and 
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j. In addition, the switchover time distribution may depend on the state of the 
last visited queue (empty or not). Remark also that, unlike in Srinivasan [15], 
we can also have JEBij{n) < W,Bij{n), which means that serving a customer 
can actually take less time than just moving. This would amount to negative 
service times in the usual polling formulation, which are not easy to handle! 

It is worth noting that our setting is valid for both discrete-time and con- 
tinuous-time models with or without batch arrivals. We do not need a separate 
analysis for each case, until the computation of waiting times, where we will 
have to give precise definitions of the vectors Bij{n) and Bij{n). We do not 
describe yet the exact arrival process and the time taken to travel from one 
station to another, since we only need the distribution of the number of clients 
arriving during such a move. However, due to the assumption of independence 
between successive arrivals, the model applies mainly to the following situations: 

• discrete-time evolution, when a batch of customers arrives at each station 
at the beginning of each time slot; batches are i.i.d. with respect to the 
time slots but they can be correlated between stations; 

• continuous time evolution, when the arrival process is a compound Poisson 
process: customers arrive at the system in batches at the epochs of a Pois- 
son process and the batches have the same properties as in the previous 

case. 

Let Xi(n) be the number of clients waiting at station i at polling instant n 
and S{n) be the corresponding position of the server. If we define X{n) = 
(Xi(n),...,Xjv(n)), then £ = {S,X) foi ins a Markov chain. Throughout 
this paper, we will assume that this chain is irreducible. For any vector z = 
{zi, . . . , zn) € (where V denotes the unit disc in the complex plane), we 
note z^ = zf^ ■ Z2^ ■ ■ ■ z^^ and define the following generating functions: 

Fiiz;n) 
F^{z;n) 

where If is as usual the indicator function of the set £. The following result 
holds: 

Theorem 1 Let A{z), A{z) and A(£) be matrices defined as follows: for all 
i,j G S, the elements of row i and column j are given by 

= Pj,iaj4z), 



= lE[z^-^(")], 

= lh[2: ^ '^{S(n)=i}\i 

= ]E[z*'"^l{s(„)=i_Xi(n)=0}]i 

= Fi{zi, Zi-i, 0, Zi+i, Zn; n), 

= lim Fi{z;n), 

= lim Fi{z;n), 

n— >-oo 
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Then the vectors F{z) = {Fi (z), ...,Fn{z}) and F{z) = {Fi{z), Fn{z)) are 
related by the functional equation 

[I - AA{z)]F{z) = [A{z) - AA{z)]F{z), (2.1) 

where AA{z) stands for A{z)A{z). 

Proof We have, for all i,jGS and n > 0, 

JE[z^(''+^H{s(n)=i,S(n+l)=j}] 

= JE[z^^"'^~^^'-^^"'^~'''llS{n)=i,Sin+l)=j,Xi{n)>0}] 

+TE[z-^^"^~^^^-^^'^H^S(n)=i,S{n+l)=j,X,{n)=0}] 
= Pi,ja^,Jiz)JE[z^'^^^^-^'t^s{n)=^,X,in)>0}] 
+Pijdij{z)lE[z^^''H[s{n)=i,X,{n)=0}] 

F,{z;n) - F^{z;n) ^ r^^r^ ^ 

= Pijo-idiz) \-Pi,jaij{z)Fi{z;n). 

Zi 

The second equality above uses the independence of the routing and the 
arrivals with respect to the past. Summing over all possible i, we get 

ix ,^Fi(z;n) — Fi(z;n) _ . . 

Fj{z;n + 1) = 2_^Pi,jaij{z) V 2_^Pi,jaij{z)Fi{z;n), 

i=i i=i 
so that, letting n — > oo, 

N N N 

p.^^ _ ^ PM^F,{z) = Y^pijaiAmi^ - E ^-'^^^m, 

1 Zi . . ^ Zi 

i—l 1=1 1=1 

which is equivalent to (2.1). ■ 

Defining, when it exists, 

D{z) = [I - AA{z)]-^[A{z) - AA{z)], (2.2) 

one sees that (2.1) can be rewritten as 

F{z) = D{z)F{z). (2.3) 

This functional equation contains all the information sufBcient to character- 
ize F{z) and F{z). Although its solution for > 3 is still an open question, 
partial results can be derived as shown in the following sections. 



3 The stationary distribution of the position of 
the server 

What renders the polling model presented in the previous section difficult to 

analyze is, among other things, the fact that the movements of the server depend 
on the state of the visited stations. In particular, {S{n)}n>o is not a Markov 
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process, as it would be when pij = pij for all i and j. Some computations are 
needed to get the stationary probabilities of the position of the server, that is 

Fde) = P{S^i,X,=Q). 

It appears that these stationary probabilities depend not only on the transi- 
tion probabilities, but also on the following mean values: 

N 



= ^Pi,j^Bij.q{n), (3.1) 

N 

= ^P-^,3^Bij-q{n). (3.2) 



Here ai-q (resp. cn-q) is the mean number of new clients at station q between 
the arrival of the server at station i which was non-empty (resp. empty) and 
its arrival at the next (arbitrary) polled station. Since the arrival process is 
stationary, it is possible to define for all g € 5 the mean number Ag of arrivals 
at station q per unit of time and write, with obvious definitions for r, and fj, 

Cti-q = XqTi, (3.3) 
0(i;q = \fi. (3.4) 

Throughout this paper we will sometimes write F (resp. F) instead of F{e) 
(resp. F{e)) in order to shorten the notation. Moreover, *F denotes the trans- 
pose of the vector F. 

Theorem 2 Define the matrices A = {oii-^q) and A = {oii;q)- When the Markov 
chain C is ergodic, F and F satisfy the following linear system of equations: 

*Fe = 1, (3.5) 
*F[I-P] = *F[P-P], (3.6) 
*F[I-A] = *F[I-A + A]. (3.7) 

Define 

N 



When p ^ 1, Equations (3.6) and (3.7) can he rewritten as 

= ^PijF'i Xiipij - Pi,j), for all j e S, (3.8) 

ies ies 

and F is given by _ 

= Fj - A,^- (3.10) 
Moreover the mean time between two consecutive visits to queue j is 

^Tjj = f (3.11) 
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The quantity f defined in (3.9) can be seen as the mean time between two 
polling instants. Note that, with the above definitions, p is in general different 
from the classical p defined in queueing theory. However, they coincide if P = P 
and the time between two polhng instants can be decomposed into a state- 
independent switchover time and a service time which only depends on the 
station where the customer is. 

Proof of Theorem 2 Although this theorem could be proved analytically us- 
ing the functional equation (2.1), we give here a simple probabilistic interpreta- 
tion of (3.6) and (3.7). Indeed, when the system is ergodic, (3.6) follows directly 
from 

JV 

P{S = j) = Yl [^(^ = i,Xi> 0)Pi,j + P{S = i,Xi = 0)pi,, 

i=l 

Moreover, writing the equality of the outgoing and ingoing flows at station 
q gives 

JV JV 

P{S = q,Xg>0) = Y^Yl = i,Xi> 0)pi,jJEBij,g{n) 

i=i i=i 

+ P{S = i,Xi = Q)pij-SEBij,g^ , 

which is equivalent to (3.7) if we take into account (3.1) and (3.2). It is straight- 
forward to see, using (3.3) and (3.4), that, when p ^ 1, 

[I - A][I - A + A]-' = I - 

1-/3 

With this relation and f as in (3.9), Equations (3.8) and (3.10) follow easily 
from (3.6) and (3.7). Let Njj-j be the number of arrivals to a queue j between 
two consecutive visits of the server. The equality of flows reads 

P{Xj>0\S = j)=-ENjj.j. 

Since the last equation can be rewritten as 

Fi - Fi 



— = — AjlETj^j , 

we obtain (3.11) and the proof of the theorem is concluded. ■ 

The set of equations (3.5), (3.8), (3.9) and (3.10) provides a convenient way 
to compute F and F. The following proposition gives simple conditions under 
which its rank is full. 

Proposition 3 The linear system given by (3.5), (3.8) and (3.9) has a unique 

solution when P has exactly one essential class. When P has K > 1 essential 
classes £i, . . . , £k, the Markov chain C is never ergodic unless 

J2 Yl ^i(Pi'j - Pi^j^ = 0' m<K. (3.12) 

je£„ ies 
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Proof Define Vi = Fi/f. Then the system (3.8)-(3.9) reads 

tpj -^Pijtjji = ^>^i{Pz,j -Pi,j), for allied, (3.13) 
ies ies 

Y^i^ih = 1-/3. (3.14) 

The system (3.13) have rank A'' — 1 unless the stochastic matrix P has several 
essential classes. Then, for each m < K, i G and j ^ £m, Pij = 0. In this 
case, (3.13) has solutions only when (3.12) holds. ■ 

The conditions (3- 12) represent in some sense "zero drift" relationships 
which, as usual, imply more involved derivations. However, when the arrivals 
form independent compound Poisson processes at each queue, the system is 
never ergodic! The proof is of analytic nature. It requires Taylor expansions of 
second order in equation (2.1), similar to those extensively used in Section 5. In 
general, when the batches are correlated, it is difficult to conclude and we sus- 
pect that all situations might occur (ergodicity, null recurrence or transience). 

It is worth noting the special role played by P. In particular, when the 
arrivals are Poisson, the Markov chain is never ergodic if P admits several 
essential classes (obviously, P has then to be chosen to ensure the irreducibility 
of jC). In our opinion, this result is not intuitive and we cannot be explained 
just by waving hands. 



4 Conditions for ergodicity 



The purpose of this section is to classify the process £, viewed as a random walk 

7N 

-+ ■ 



on >S X Z^, in terms of ergodicity and non-ergodicity. 



4.1 Necessary condition 

We give a necessary condition for ergodicity which is a simple consequence of 
the results of Theorem 2. When the system is ergodic, we have Fi > for all 
i G S and Equation (3.10) implies the following. 

Theorem 4 // the Markov chain C is ergodic, then 

p<l, (4.1) 

Ajf < Fi, for all i G S. (4.2) 

These conditions extend in the state-dependent case the results obtained 
in [1]. When they hold, we prove another useful inequality. Indeed, instantiat- 
ing (4.2) in (3.9) yields 

1-p 

or, equivalently, 

1-EV.>0. (4.3) 
jes 
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4.2 Sufficient condition 



Let lis emphasize that this part could be skipped by readers not strongly in- 
terested (if any!) in ergodicity. The mathematical understanding requires the 
reading of Appendix A.l, which summarizes deep works of [12], [4] and [5]. The 
approach relics on the study of a dynamical system associated to £. For the 
sake of readability, we recall hereafter two main notions. 

• Induced chain . It is a Markov chain corresponding to a polling system 
in which the queues belonging to the face A C S are kept saturated. From 

a purely notational point of view, the original system would correspond 
to the case A = 0. The behaviours of £ and £^ are not directly connected 
to each other. 

• Second vector field v^. One can imagine that the random walk starts 

from a point which is close to A, but sufficiently far from all other faces 
A', with A (Z! A'. After some time — sufficiently long, but less than the 
minimal distance from A' — , the stationary regime in the induced chain 
will be installed. In this regime, one can ask about the mean along A: it 

is defined exactly by v^. 

As in Theorem 2, the stationary position of the server for any ergodic induced 
chain can be obtained as the solution of a linear system consisting of A'' + 1 
equations. Indeed, for all t G S, j ^ A, we have 

n^it) = Y.n^{s)ps,t + Y.^^i^)Ps,t (4.4) 

seA s^A 

+ ^Tr^{s,Xs > 0){Ps,t -Ps,t), 
Tr^{j,Xj>0) = 5^7r^(s)AjT, + 5^7r^(s)Ajf, (4.5) 

sGA s^A 

+ ^TT^{s,Xs > 0)Aj(Ts -fs), 
S0A 

J2n^is) = 1. (4.6) 



Let us introduce the quantities 



sgA 
1 



^7r^(s)T,+^7r^(s)f, 

'-seA S0A 

Note that when A is ergodic, -K^{j,Xj > 0) = Ajf^, for all j ^ A, and 
f ^ can be interpreted as the mean time between two polling instants of the 
induced chain. Then the system defined by (4.4) and (4.5) can be replaced by 
the forthcoming A'' + 1 equations: for all t € S, 

n^it) = 5^7r^(s)p,,t + ^7r^(s)p«,t + f^5^A,(p,,t-p,,t), (4.7) 

sGA s^A s^A 

1 



5^7r^(s)r, + ^7r^(s)f, 

■s6A s^A 



(4.8) 
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When A = 0, the system (4.6), (4.7) and (4.8) coincides with (3.5), (3.8) 
and (3.9). In addition, under (4.3), we have 



< 1, for all A 's. 



(4.9) 



An easy flow computation shows that the components of the drift vector 
M{s,x) can be expressed as 

Mj{s,x) = XjTsl{x^>0} + ^jfs'^{xs=0} - '^{xs>0,s=j}- 

Then the computation of the second vector field becomes easy. For j £ A, 

xec^ 
ses 

= ^ ■JT^is,x)Mj{s,x) + ^ TT^{s,x)Mj{s,x) 



s0/\ 



sGA 



= XjY^\TT^{s,Xs><d)Ts+TT''{s,Xs^<d)i 
sgA 

+ A,-^7r^(5)T,-7r^(i) 

= A,f^-7r^(j). 
In order to apply Theorem 11, it will be convenient to introduce 



fiix) 



'x[I-A + A]-^ 



N 



Xi + 



(4.10) 



This function is not as outlandish as it could seem at first sight. It is indeed 

directly related to flow conservation equations (see [5]). Let A be any ergodic 
face. For i € A, we get after a straightforward computation 



A,: 



N 



N 



1-/5 



For i ^ A, and from the very definition of fi{x), we have 



N 



or 
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Theorem 5 Assume that (4-l)-(4-2) hold and that, for any ergodic face A, 

N 

fi{v^) = ^^i^^i^^^J^ -Tr'^d) <0, for all i€ A. (4.11) 

Then the random walk C is ergodic. In particular, when P = P , the condi- 
tions (4-1) and (4-2) are necessary and sufficient for the random walk £ to be 
ergodic. 

Proof Wc shall apply Theorem 11 (quoted in Appendix A.l), which in itself 
contains a principle of gluing Lyapounov functions together. Most of the time, 
these functions are piecewise linear. Since the fi{v^) enjoy "nice" properties, it 
is not unnatural to search for a linear combination like 

JV 



f{x) =J2uifi{x), 

i=l 



where u = {m, . . . ,ujv) is a positive vector to be properly determined. Then, 
for any ergodic A, 

/(.-) = J:[u, + %^(r. - f.)]/.(.-). 

iGA ^ 

The basic constraint for u is to ensure the positivity of /. Using (4.3) 
and (4.9), it appears that a suitable choice is Ui = max(fj — Ti,e), for some 
e, positive and sufficiently small. Then one can directly check that f{x) > 0, 
for any x G K^, and f{v^) < 0, for any ergodic face A. ■ 

Although the conditions of Theorem 5 are sufficient for crgodicity, they might 
well be implied by conditions (4.1)-(4.2) only. We did not check this fact, since 
it is difficult to compare algebraically the solutions of system (4.7)-(4.8), for 
different A's. Let us simply formulate the conjecture that (4.1)-(4.2) are also 
sufficient for crgodicity when P ^ P. 

Theorem 6 Assume that P = P and that one of the following hold: 

(i) there exists j € S, such that Xjf — Fj > 0; 

(m) p > 1. 

Then the random walk L is transient. 

Proof We assume that the queues are numbered according to the following 
order: 

Consider all faces Aj = {i, . . . , A''}, for i = 1, . . . , A^. The idea is to show 
the transience of C by visiting the ordered set Ai, A2, . . . Aat. In the usual 
terminology of dynamical systems, this amounts to finding a set of trajectories 
going to infinity with positive probability. To this end, the following algebraic 
relationship will be useful: for any A and k ^ A, we have 

(Afcf^+W - Ffe)(l - p'^+{'=>) = (Afef^ - Ffe)(l - p^). (4.13) 
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By definition, tlic face Ai = 5 is crgodic (sec Appendix A.l). Assume now 
Aj is ergodic, for some fixed i £ S. If v^^ < 0, tfien one can prove tlie inequality 



which in turn yields 1 — p^' > 0. By using (4.13), (4.12) and Theorem 5, it 
follows that the face Aj+i is ergodic. Conversely, if v^' > 0, then v^^ > and 
the random walk C is transient (see [4]). 

Thus, by induction, we have shown that cither C is transient or the face 
Ajv = {N} is ergodic. In the latter case, the assumptions of the theorem, 
together with (4.13) and (4.12), yield v%'' = Xnt^" - Fn > 0, so that £ is 
transient. ■ 

4.3 Remarks and problems 

(i) In general, the second vector field requires to compute invariant mea- 
sures of Markov chains in dimensions strictly smaller than N (the induced 

chains). This is rarely possible for TV > 3, but there are some miracles, 
the most noticeable ones concerning Jackson networks and some polling 
systems [5]. 

(ii) No stochastic monotonicity argument seems to work, when the Tfc — ffc's 
arc not all of the same sign, even for P = P. This monotonicity exists 
and is crucial in [1], [7] and [8]. 

(ill) The analysis of the vector field is interesting in itself and will be the 
subject of a future work. Let us just quote one negative property: when 
A is ergodic and < 0, then Ai = A \ {k} can be non-ergodic, contrary 
to what happens in Jackson networks (sec [5]). In the case P = P, we 
conjecture nonetheless that the dynamical system formed from the velocity 
vectors is strongly acyclic (in the sense that the same face is not visited 
twice). 

5 The first moment of the queue length 

This section shows that, when the system enjoys some symmetry properties, it 
becomes possible to derive the stationary mean queue length ]E[Xm \ S = m] 
seen by the server at polling instants. To this end, ergodicity of £ will be 
assumed, as well as the existence of second moments of all random variables of 

interest. 

Assumption Ai The system has a rotational symmetry 

(i) for all i,j € S, ptj = pj-i, pij = pj-i, aij{z) = aj-i{z) and dij{z) = 
aj-i{z); 

(ii) for ain<s,d< N, 




Udizi, . . . ,zn) = aci{zi+. 
ad{zi, . . . ,zn) = ad{zi+, 



. . . , Zn^s), 
■ ■ ■ , Zn+s)- 
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Note that this assumption is less restrictive than the ones appearing in the 
literature, even for state-independent situations. It allows in particular cyclic 
and random polling strategies. From Ai, we directly get, for all 1 < i, s < N, 

Fi{z) = Fi+s{zi+s, - ■ ■ ,ZN+a), (5.1) 

Fi{z) = Fi+s{zi+s,---,ZN+s)- (5.2) 

Consequently, P{S = i) = Fi{e) = 1/N. In addition, A{z) and A{z} are 
circulant matrices — see Appendix A. 2 for the results and notation which will 
be used from now on. It implies in particular that, for k G S,Vk is an eigenvector 
of A{z) and A{z} with respective eigenvalues 

N 

l^k{z) = ^Pdad{z)u!'^, 

d=l 

N 

l^k{z) = ^Pda,d{z)uj'l. 

d=l 

Note that ij.n{z) (resp. fLN{z)) is the generating function of the number of 

arrivals during the transportation of one client (resp. a move without client). 

Moreover, /x^ = iJ-k{e) and jlk = /ift(e) are eigenvalues of P and P. For a priori 
technical reasons we shall also need 



Assumption A2 For any k < N , flk ^ I- 

These relations are not surprising in view of Proposition 3. Since the system 
is symmetrical, the mean number of customers arriving during a travel of the 
taxicab ai-q and cti-q as defined by (3.1)-(3.2) does not depend on i and q; we 
note them a and a and remark that 

The second derivatives of ij,n{z) and P,n{z) are defined as 

~ dZqdZr ^ ^ ' 

(,) def d^flNi^ , . 

"'^ dZqdZr ^ ' ^ ' 

Note that these values depend only of the unsigned distance between q and 
r and that we can write a^^'^^^^ = af_^^^^ = af, a^^'^_^^ = a^![^ ^ = af. It will 
be also convenient to introduce tlic first and second moments of the number of 
arrivals between two polling instants, that is 

a = P{Xm >0\S = m)a + P(X„ = | 5 = m)a, (5.5) 
a'^' P(X„>0|S = m)af +P(X„ = 0|5 = m)af . (5.6) 
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Using the symmetry of the system and according to the properties of gener- 
ating functions, we have for any m € S 



8F 



(5.7) 



By symmetry, this expression is independent from m. In order to derive 
]E[Xj„ I 5* = to], we use the fact that the matrices / — A/S.{z) and A{z) — AlS.{z) 
are not invertible when z = e. This imphes that the matrix D[z) as defined 
in (2.2) is not continuous in a neighborhood of ^ = e and provides a mean to 
compute the derivatives of F{z) and F{z) for z — e. However in our setting, 
even the existence of D{z) in a neighborhood of e is difficult to prove and the 
computations become really involved. The approach that we present here avoids 
the theoretical problems and simplifies the computations. In the next lemma, we 
show how it is possible to choose z such that the left member of Equation (2.1) 
can be rendered of "small" order w.r.t. the right one. 

Lemma 7 Let t ^ Zt be a function from IR_ — ?■ such that in a neighborhood 
oft = 0, we can write Zt = e + tz + \t^z + o(t^), where z = z\e\-\ h ijvejv = 



Ci"?! + • • --^CnVn € C"'^ and z = (zi, . . . , z^r) G C"'^. Assume that zi + • 
0. Then there exists for t > a vector Ut such that, for 1 < k < N , 



■ + zn 



*ut [I - AA{zt)\vk = 
*ut [A{zt) - AA{zt)]vk = 



o{t), 



CiV-fe + o{t), 



and 



*ut [I - AA{zt)]e 



N-l 



CiC 



N-l 



1 N ^ 



(5.8) 
(5.9) 



,1-Na 
2N 



N 



^Zg + 0{t''), 



(5.10) 



9=1 



*ut [A{zt) - AA{zt)]e 



1=1 
,1 - 



JV-; 



N 



+ ^y-(a<^'-a<=')a;f 



2 

d=l 



Na + Na .. , 



IN 



(5.11) 



9=1 



Proof See Appendix A. 3. 



For the sake of simplicity, the computations have been carried out in the 
most natural way which apparently requires /ifc ^ 1 for fc < N . In fact, the 
reader can convince himself that all derivations could be achieved by rendering 
the right member of (2.1) small w.r.t. the left member. This would yield the 
same result without any restriction on /Xfc, at the expense of a somewhat longer 
proof. 

One problem in polling models with 1-limited service strategy is that the 
rank of the systems of equations that we can write is [n/2j, which means that 

dF / dzm{e) cannot be computed. Fortunately, under the following assumption, 
we are able to compute the most important value, that is dFm/dzm{e)- 
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Assumption A3 The routing matrices are such that 



[I-P][I-*P] = [I -*P][I - P]. 

It is not easy to describe all models satisfying A3. The simplest one is 
the classical Markov polling obtained when P = P, for which the results of 
Theorem 8 thereafter are greatly simplified. Another admissible model is when 
the routing probabilities depend only on the absolute distance between stations, 
that is when P = *P and P = *P . 

Theorem 8 Assume that A\, A2 and A^ hold. Then, for any station m, the 
stationary probability that the server finds station m empty is given by 

P{Xm = 0\S = m)= . . (5.12) 

1 — Na + Na 

Moreover the mean number of customers found by the server when it arrives 
at station m is obtained by means of the following formula 

Naa ^ ^ 



nXm\S = m] = Na+-^J2T^ 

1 - Na ^ 1 - 

1-Na + NaN Na-Nal,^ 
+ 1-Na T"^ + 1-Na 2^"'' 



N-l . N 

1-Na^^ 1-W2^i ' ' 

Finally, the mean number of customers present at arbitrary station m at a 
polling instant can be expressed as 

N-l 

JEXm = a+- ^ 

1 - Na 1 - A*; 

N 

l-Aa + AaA_(2) Na - Na I ^-^ 
+ 1-Na y + 1-Na 2^"'^ 

N — 1 N 

Proof See Appendix A. 3. ■ 



6 The mean waiting time of a customer 

Theorem 8 is the key to compute of the waiting time of an arbitrary customer. 
In this section, we focus on the continuous-time case with compound Poisson 
arrivals and further assume that the arrivals at each station arc independent. 
While Theorem 8 is valid in discrete time or with correlated arrivals, we do not 
give any formula in these cases, since the notation would become too involved. 

The notation is as follows: customers arrive in i.i.d. batches of mean length 
b and second moment 6*^' at the instants of a Poisson stream with intensity A 
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at each station. The intensity of the arrival process at each queue is A = A6. 
The time between two polling instants when the first queue is not empty (resp. 
empty) has mean and second moment r and r'^' (resp. f and f'^'). We must 
keep in mind that in general t depends on P like ai-q in Equation (3.1): if 
Td is the mean time between two consecutive polling instants at stations i and 
i + d, then we have r = pm + • • • + PnTn- However, in most symmetrical 
polling models that have been previously analyzed, the switchover times are 
equal and r does not depend on P. We will say in this case that the stations 
are equidistant, which is obviously not the case in our original taxicab problem. 
The same holds for r'^' , f and f'^' . 

Theorem 9 Assume that Ai, A2 and As hold. Then the mean stationary wait- 
ing time of a customer is given by the relation 



JE[W] 



1-NXt 



N-l 

E 



+ ■ 



6'=' - b 



26 



1 



1 - W 



NXr- 



2{1 -NXt) 2f 



T+{N-l)f 

1-N\t 



1 - TV At ^ 



(6.1) 



Proof See Appendix A. 3. 



One surprising consequence of this formula is that when the arrival process is 
Poisson and the stations are equidistant, the mean waiting time does not depend 
on the routing P. In other words, it does not depend on where the server goes 
after serving a customer. 

Beside giving the mean waiting time for a customer in our taxicab model, 
Equation (6.1) contains in fact many known formulas for waiting times corre- 
sponding to various service and polling strategies. In the next subsections, we 
give some of these applications. 

A notation closer to the classical qiicuing theory notation is necessary to 
make the link with known results. Let w (resp. w) be the mean walking time — 
or switchover time — from a non-empty (resp. empty) station and let a be the 
mean service time required by the customers. We denote by w^^^ , ti)'^' and cr*^^' 
the associated second moments. 



6.1 A state-independent Markovian polling model 

The first possible application of our model is the classical state independent 
polling model with 1-limited service strategy. In this model, we have P = P, 
T = w + a, f = w and (6.1) becomes 



NX{w + a) f-^ I- m 2{l-NX{w + a)) 2w 
6<=^> - 6 Nw + a 



2b 1 - NX{w + a) ■ 



This formula is valid for all symmetric polling models with Markovian polling. 
It is interesting to note that, in the case of equidistant stations, there is only 
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one term of the formula which depends on the routing. To compute this term, 
we remark that, if P{x) is the characteristic polynomial of P, 



N-l ^ N-1 



-1 - w 



^ 1 

— log 
dx 



V{x) 



x-1 



x=l 



The case of the cyclic polling is obtained by taking pi = I and Pd = for 
1. Then the eigenvalues of P are oJi, . . . ,ujn and P{x) = x^ — 1. So in this 
case, we have 

JV-l 



1=1 



_ N-l 
•W ~ 2 



Another classical polling strategy is the random polling with = 1/N for 
all d. In this case, we have m = ioi I < N and 



= N -1. 



The comparison between formulas for cyclic polling and random polling 

shows that the mean waiting time is smaller in the cyclic case. In fact this 
property can be generalized to all Markovian polling models. 

Lemma 10 Among all Markovian polling strategies for symmetrie 1-limited 
polling models with equidistant queues, the cyclic polling strategy minimizes the 
waiting time of the customers. 

Proof This result is very easy to prove from Equation (6.1). We have 

JV-l , JV-l 



2(i-K(mO) + Im;P-i' 

where ^{z) is the real part of z. Since < 1 and < 1, we find that 

JV-l 



1=1 



N - 1 
— > . 



The bound is attained if and only if for alH, = 1. Since is the center 
of gravity of a; ... , cj/^ with weights pi, . . . ,pi\[, this is only possible when for 
some s G <S we have Ps = 1 and = for d ^ s. When s is not a divider of 
A'' (to ensure that the routing matrix is irreducible), this is equivalent to cyclic 
polling. ■ 

This property also holds for discrete time systems and systems with corre- 
lated arrivals. The computation would be more difficult in the case where the 
stations are not equidistant. 
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6.2 The exhaustive and BernouUi service strategies 



An iiitcrcisting application of our model is to show that sonic service strategics 
can be obtained by a proper choice of polling strategy. In this subsection, we 
show how we can apply our model to exhaustive and Bernoulli service strategies. 

Consider the case where the stations are equidistant, except that the switchover 
time from a station to itself after a service is zero. Then, if we denote pn = 1 — tt, 
we have r = nw + a, r'^^ = nw'-^^ + 2-kwcf + a^^^ , f = w and f'^' = w'-^^ . If we 
assume for the sake of simplicity that the arrival process is Poisson, (6.1) yields 



JV-l . 

^ ^ 1 - NXinw + u) 4-1 l-n 



w 1 7VA(7rw'^' +27r«;<T + CT("') 



w 



(2) 



l-NX{-KW + a) j^^l-jli^ 2{l- N\{-KW + (j)) ^ 2w ' 

The value of IE[iy] depends on P only through the value of p^r. One known 
model that is described by this equation is the Bernoulli strategy of Servi [14]: 
when the server has served a customer, it quits the queue with probability tt and 
continues to serve it with probability 1 — tt. In the original model, the server 
polls the queues in cyclic order. Here, we have a Markovian Bernoulli polling 
model if we take P = (1 — tt)/ + ttP. It is easy to check that this choice of P 
satisfies A3. Moreover, IE[H^] is an increasing function of tt that is minimal when 
TT = 0. TJiis case corresponds to an exhaustive service strategy with Markovian 
routing P: the server polls the same queiie until it is empty and then moves to 
another queue using the routing matrix P. In this case. Equation (6.1) simply 
becomes 

^ ^ ~ 1 - NXa ^ 1 - w 2(1 - NXa) ^ 

Using Lemma 10, we find that the above expression is minimal when P 
describes a cyclic polling scheme. Further comparisons between different polling 
strategies can be found in Levy, Sidi and Boxma [11]. 



Appendix 

A.l Second vector field 

This appendix contains some basic definitions and results from the theory of 
dynamical systems. These definitions have been adapted for the Markov chain 
C defined in Section 2. We refer the reader to [12] and [4] for a more complete 
treatment of the subject. 

Faces. For any A C {1, ... , N), define C 1R+ as 

B'^ = {{xi,...,Xn) : Xi > 0<(^i G A}. 

-B^ is the face of associated to A. Whenever not ambiguous, the face 
B'^ will be identified with A. 

Induced chains. For any A 7^ <S we choose an arbitrary point a e B'^ n 2^ 

and draw a plane of dimension iV — | A | , perpendicular to B^ and containing 
a. We define the induced Markov chain with state space 5x (C^n2^) (by 
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an obvious abuse in the notation, we shall write most of the time S x C^) and 
transition probabihties 

AP{s,x){t,y) =P{s,x){t,y} + ^^P{s,x){t,z), , ^ S , S,t G S, 

where the siimmation is performed over all z G 2^ , siieh that the straight line 
connecting z and y is perpendicular to . It is important to note that this 
construction does not depend on a. 

Assumption A4 For any A the chain is irreducible and aperiodic. A is 
called crgodic (resp. non-ergodic, transient J according as jC^ is ergodic (resp. 
non-ergodic, transient). 

For an ergodic let 7r^(s,£'), (s,x) £ 5 x be its stationary transi- 
tion probabilities. Moreover, let 7r^(s) (resp. it^{s,Xs > 0)) be the stationary 
probability that the server is at station s (resp. and finds it nonempty). 

First vector field. For any (s, x) ^ S x 'Z.^^ the first vector field is simply the 
mean drift of the random walk C at point (s, x): 

M{s,x)= ^ {y-x)p\x{n+l) = y,S{n + l) = r X{n) = x, S{n) = s 

res, ye^^ 

Let = {vi, ... ,v'^) such that 

Vi = 0, i ^ A, 

= ^ 7r''(s,f)M,(s,f), i G A. 
(s,x)eSxCA 

Intuitively, one can imagine that the random walk starts from a point which 
is close to A, but sufficiently far from all other faces B'^ , with A ^ A'. After 
some time (sufficiently long, but less than the minimal distance from the above 
mentioned B'^ ), the stationary regime in the induced chain will be installed. In 
this regime, one can ask about the mean drift along A: it is defined exactly by 
. For A = iS, we call A ergodic, by definition, and put 



sG5 



From now on, when speaking about the components of ■u^, we mean the 
components with i G A. 

Assumption A5 7^ 0, for each i € A. 

Ingoing, outgoing and neutral faces. Let us fix A, Ai, so that A D Ai, A 7^ 
Ai, that is to say B D B^^ {B is the closure of B^). Let B^ be ergodic. 
Thus is well defined. There are three possibilities for the direction of 
w.r.t. B^'^ . We say that B^ is an ingoing (resp. outgoing) face for B^^ , if all 
the coordinates for i G A\ Ai are negative (resp. positive). Otherwise we say 
that B^ is neutral. As an example we give simple sufficient criteria for a face 
to be ergodic. 
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The second vector field. To any point x € IR_|_ , we assign a vector v{x) 
and call this function the second vector field. It can be multivalued on some 
non-ergodic faces. We put, for ergodic faces B^, 

If B^^ is non-ergodic, then at any point x G B^'^, v{x) takes all values 
for which B^ is an outgoing face with respect to B^^ . In other words, for x 
belonging to non-ergodic faces, with ||^|| sufficiently large, 

x + v{x) e IR+ , 

for any value v{x). If there is no such vector, we put v{x) = 0, for x G B^'^ . 
Points X e R^, where v{x) is more than one- valued, are called branch points. 

There are few interesting examples, for which only the first vector field suf- 
fices to obtain ergodicity conditions for the random walk of interest, but it is 
nevertheless the case for Jackson networks. In general, the second vector field 
must be introduced as shown in the following theorem, taken from [12] (and 
extended in [5] to the case of upward unbounded jumps). 

Theorem 11 Assume that there exists a nonnegative function f : IR^ — >■ 1R+ 
such that 

(i) for some constant c > 

f{x)-f{y)<c\\x-y\\; 

(m) there exists 5 > 0, p > such that for all x G B^ 

f{S+v{x))- f{x)<-8. 

Then the Markov chain C is ergodic. 
A. 2 Some simple results about circulant matrices 

In this appendix, we recall some well-known properties of the circulant matrices 
which are used throughout this paper. These properties are given without proof, 
since they can easily be verified at hand. Throughout the paper, we take the 
convention that every subscript less than 1 or greater than N should be shifted 
into the correct range. 

Let (ei, . . . , e/v) denote the canonical basis of and let e = (1, . . . , 1). 
Moreover, we define ujk == exp(2t7rfc/A^), with = —1. 

Definition A circulant matrix M of size N is a matrix whose coefficients rriij 
verify the relation: 

mi+k,j+k = for all 1 < jj ^ < 

The following Lemma summarizes some key properties of these matrices: 



19 



Lemma 12 Let M be a circulant matrix of the form: 



mi mjv 



mi \ 



M = 



\ ruN-i mN-2 



ruN j 



Then for each 1 < k < N, the vector Vk 



= X^i^i w^'e, is an eigenvec- 



N 

tor of the matrix M with eigenvalue X^^^i u\mi. Moreover, {vi,. . . ,vn) is an 
orthogonal basis of in which ei, . . . , ejv can be expressed as 



One important feature of circulant matrices is that they share the same basis 
of eigenvectors. Note that, with the notation given before Lemma 12, we have 
vn = e. Circulant matrices enjoy other properties that are not used here: for 
example, the product of two circulant matrices is a circulant matrix and this 
product commutes. 

A. 3 Proofs of the results of Sections 5 and 6 

The derivations in this section are essentially analytic. 

Proof of Lemma 7 This proof uses specific properties of circulant matrices 
to perform a fine analysis of the behavior of A{z), A{z) and A(z) in the neigh- 
borhood of i* = e. As pointed out later in the proof, we study these matrices 
only for z G , thus avoiding any analytical continuation. The basic relation 
used thereafter is a simple consequence of the definitions of A{z) and A(z): 



where A.q{z) stands for the q-th column of A{z) and can be written as 




fe=i 



1 

AA{z)vk = Hk{z)vk + ^( l)ujl'^A.q{z), 



(A.1) 




i=l 




N 




So, if we define 
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(A.l) can be rewritten as 



N 



[I-AA{z)]vk = (^-2) 

1=1 

N 

[A{z) - AA{z)]vk = {fik{z} - Hk{zj)vk -^£i-k{z)i^i{z)vi (A.3) 



Since *e[I - AA{e)] = *e[A{e) - AA{e)] = 0, we define, for any set of 
arbitrary complex numbers (ci, . . . , cjv-i), the vector Ut as follows: 

1 / 

1=1 

With this definition, we have 

*- , N-l 

*ut [I - AA{z)]e = -^[I - AA{z)]e+ - ^ q V; [/ - AA{z)]e 



1=1 



N-l 



1 - /Xiv(z) - eo{z)idN{z) - t ^ ci£i{z)iii{z). (A.4) 



When t is small, we have the relations: 



N ^2 N 

HN{zt) = l+taY^Zq+ — a^q,rZqZr 
q=l ^5'^=! 



J.2 



2 ,=1 

AT ^2 ^ 

/iAr(f() = l+ta^Zg + — ^ a'g,l-ZqZr 
q=l <7)'^=1 



J.2 ~ -'^ 



2 «=1 

AT ,2 Af 



q=l q=l 
9=1 

The formulas given in Appendix A. 2 allow to relate ii, iiv to Ci, ... , C,n'- 



g=l 1=1 q=l 1=1 

Using the remark after the definition of a^^^ and a^^^, we find that: 

<r^qir = E E = E c^Civ-; E ^>f- 

q,r=l d=l q=l 1=1 d=l 
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Applying (A. 5) and (A. 7) to (A. 4), wc sec that the first-order term in the 

expression (A. 4) is {a — 1/N)t[zi + h ijv]. So a necessary condition to have 

a formula like (5.10) is ii H h ijv = 0. Using this relation, we have 



*Ut[I-AA{zt)]e = i'^CiOw-t'^OC 



JV-l JV 



1=1 1=1 

1=1 d=l 

+ t''-^f:z, + o{t'). (A.8) 

9=1 

We have to check that it is possible to have zt £ and zi + ■ ■ ■ + zpf = 0. 
One easy way to satisfy these constraints is to ensure that all coordinates of Zt 
arrive to 1 tangentially to the unit circle as t goes to 0. This is the case when, 
for any q, Zq is an imaginary number and Zq < 0. Moreover, for k < N, 

*Ut[I - AA{zt)]Vk = -SN-k{Zt)lJ.N{Zt)+tCk{l- lJ,k{zt)) 

N-1 

ciei-k{zt)iJi{zt) 

1=1 

= tCN-k+tCk{l-fik)+o{t). (A.9) 

This shows that we get equations (5.8) and (5.10) from (A.8) and (A.9) if 
we take 

CN-k 

Ck = -z • 

1 - Mfc 

With this choice of Ci, . . . , Cjv-i, we obtain equations (5.11) and (5.9) in 
exactly the same way. ■ 

Proof of Theorem 8 The idea of this proof is to apply the results of Lemma 7 
to the equation 

*ut [I - AA{zt)]F{zt) = *ut [A{zt) - AA{zt)]F{zt). (A.IO) 
Define Zt as in Lemma 7. Then 

N 

dZr 



BF 

F{zt) = F{e)+tY,ir^{^ + o{t)- 



Moreover, Equation (5.1) implies that dFi/dzr{e) = dFi-r+m/dzm{e) for 
any fixed m £ S and 

OF ^9F 



i=l "^"^ -'^ k=l 
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N 



= E 



k=l 



''^ dZr, 



Wj, - — [e)Vk, 



where {ipi{z), . . . , (^^{z)), defined as 



N 



i=l 

are the coordinates of F{z) in the basis {vi,- ■ . ,vn). Finally, 



F{zt) = F{e)+tNj2Ck^rQ^i^Vk + 5{t). (A.ll) 



k=l 



With a similar definition for <^{z), 



N 



fe=l 



F{zt) = F{e) +tNY, Cfe^fc '"l^leliTfe + 3{t). (A.12) 



We now apply Lemma 7 and Equations (A.ll) and (A.12) to Equation (A. 10) 
and use the fact that, by symmetry, F{e) = Fi(e)e and F{e) = Fi{e)e 



( ^-^ r 1 N ^ 



(2), ,d 



Mi 



1 - Na 



N 



2N ^ « 

9=1 

N 



I ;=i ^' ^ d=i 



+ — 2^^ 



2^^ 

Af-l 

+ t^N E CfeCiv-fc^r 

/c=l 



g=l J 

1 - /ffc (j^fc 
1 - /ifc 



(e} + o(i2). 



As z can be chosen freely, we can derive a first equality from this equation, 
namely 

l-ATa I- Na + Na 



2iV 



2iV 



Since P(X„ = | 5 = m) = Fi(e)/Fi(e), this yields Equation (5.12). 
Taking in account (5.6), we find 



iV-l 



Ell „l-Mfe,.5</?fe 
CfeCiV-feW^ — A^:— (e) 

fe=i 



JV-l 

= - E CiCiV-i 



AT 



Fi(e)-i^i(e) ^ d 



1 - M( 



(A.13) 
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This equation contains in fact a system of linear equations that can be built 
by choosing i. The problem is that in general the rank of this system is only 
[n/2j. However, since fik (resp. fJbN-k-, P-k, pN-k) is the eigenvalue of *P (resp. 
P, *P , P) associated to the eigenvector Vk, A3 implies that, for 1 < fc < A'', 



e IR. 



As noted in the proof of Lemma 7, zi,. . . ,Zn are imaginary numbers and 
CfcCiv-fc is a negative real number for any k < N. Hence we can choose in (A. 13) 



CfcCiV-fc = -min 



The combination of the resulting equation with a similar equation containing 
only the terms where (1 — /ife)/(l — /ife)<0 yields 



N-l 



fe=i 



dVk 

dZm 



JV-l 

El 



Fi(e)-i^i(e) , 1 



AT 



1 - w 



-(2) . 



(A.14) 



We know that, by definition, Fm{z) does not depend on z^rf Hence, 

N 



dZr, 

and, using Equation (A.14) 

N-l 



(e-)=E 



uJu - — (e) = 0, 



fe=i 



Fi(e)-i^i(eO , 1-M;1 



l-fii2 



d=l 



(A.15) 



The last step of the demonstration is to get dF„i/dzm(e) from these results. 
This can be done as in Lemma 7, but with ut = e/N and Zt chosen differently. 
Using Equations (A. 2), (A. 3) and (A.5)-(A.7) with the same notation as in 
Lemma 7, we have for k < N 



N 



[I - AA{zt)]e = 1 - i^Nizt) - eo{zt)iiN{zt) 



N 



^[A{zt)-AA{zt)]e 



N 



t{l - Na)CN -t^Y. CiCiv-; + i'A^aC; 



N 



1=1 

N 



1=1 



[I-AA{zt)]vk = 



-eN-k{Zt)flN{zt) = tCN-k + o{t) 

m{zt) - m{zt) - eo{zt)nN{zt) 
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N 



t{l -Na + Na)CN -t'Y. + ^'^"C^ 



1=1 



,N 



N 



N 



1=1 



N 



[A{zt) - AA{zt)]Vk = -SN-k{Zt)lJ.N{Zt) = KN-k + 0{t) 



Combining these equations with (A. 11) and (A. 12), we find 



AT 



N 



>2) d 

d 



. 1=1 

N-1 



d=l 



NaC 



N 



k=l 



N 



N 



. 1=1 



NaC% 



N-l 



k=l 



-(e) 



+ t^N{l -Na + Na)C%^{e) + o{t^). 



(A.16) 



This equation in turn gives for (1 — Na)^^ = 1 and QCn-i = li{l<l<N 



N 



k=l 



1 - Na dz, 

N-l N 



N-l N N 



1=1 d=l 



d=l 



This expression, together with (A. 15) and (5.7), yields Equation (5.13). The 
derivation of (5.14) is done in a similar way: we use the fact that 



N 



We use Equation (A.16) with (1 - Na)C% = 1 and CfeCjv-fe = for fc < A'' 



and find 
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This gives Equation (5.14) and concludes the proof of the theorem. 



Proof of Theorem 9 We see easily that the mean numbers of clients arriving 
between polling times are respectively a = At and a = Xf. Moreover, using 
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Equations (5.3) and (5.4) and the properties of generating functions we find 
that, for d < iV 



(2) 



2^(2) 



= a'=' + A(6<"' - 6)r 



4' - 



^,(2) _ -(2) 



'AT 



a'' 



+ A(6<^' -&)f. 



The computation of waiting times uses the following classical argument: a 
non empty queue visited by the server can be decomposed into 

• the head of line customer; 

• the clients who arrived after him during his waiting time; 

• clients who arrived in the same batch as the first client, but are not yet 
served; since, by renewal arguments, the mean size of this batch is jh, 
the mean number of remaining clients is (6'^' — h)/2h. 



This can be written as 



]E[X„ I 5 = m,X„ > 0] = 1 + X^[W] + 



6<"> - 6 
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and, using the fact that X^^^i '^d^i — '^n ~ 
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This gives Equation (6.1). 
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